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AN  OPTIMAL  ALGORITHM  FOR  FINDING  THE  KERNEL  OF  A POLYGON 


D.  T.  Lee  and  F.  P.  Preparata 
University  of  Illinois  at  Urbana -Champaign 

Abstract 

The  kernel  K(P)  of  a simple  polygon  P with  n vertices  Is  the  locus  of 
the  points  Internal  to  P from  which  all  vertices  of  P are  visible. 
Equivalently,  K(P)  Is  the  Intersection  of  appropriate  half-planes  determined 
by  the  polygon's  edges.  Although  the  Intersection  of  n generic  half- 
planes Is  known  to  require  time  O(nlogn),  we  show  that  one  can  exploit 
the  ordering  of  the  half-planes  corresponding  to  the  sequence  of  the 
polygon's  edges  to  obtain  a kernel  finding  algorithm  which  runs  In  time 
0(n)  and  Is  therefore  optimal. 


This  work  was  supported  by  the  National  Science  Foundation  under  Grant 
NSF  j4CS  76-17321  and  by  the  Joint  Services  Electronics  Program  under 
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AN  OPTIMAL  ALGORITHM  FOR  FINDING  THE  KERNEL  OF  A POLYGON* 

D.  T.  Lee  and  F.  P.  Preparaca 

Coordinated  Science  Laboratory 
University  of  Illinois 
Urbana,  Illinois  61801 

1.  The  kernel  K(P)  of  a simple  polygon  P Is  the  locus  of  the  points 
Internal  to  p which  can  be  joined  to  every  vertex  of  P by  a segment  totally 
contained  In  P.  Equivalently,  If  one  considers  the  boundary  of  P as  a 
counterclockwise  directed  cycle,  the  kernel  of  P Is  the  Intersection  of  all 
the  half-planes  lying  to  the  left  of  the  polygon's  edges. 

Shamos  and  Hoey  [l]  have  presented  an  algorithm  for  finding  the  kernel 
of  an  n-edge  polygon  In  time  O(nlogn) . Their  algorithm  is  based  on  the  fact 
that  the  Intersection  of  n generic  half-planes  can  be  found  in  time  O(nlogn) ; 
they  also  show  that  O(nlogn)  Is  a lower-bound  to  the  time  for  finding  the 
Intersection  of  n half-planes.  However,  this  lower-bound  does  not  apply  to 
the  problem  of  finding  the  kernel,  since  In  the  latter  case  the  half-planes 
are  ordered  according  to  the  sequence  of  the  edges  of  P,  nor  does  their 
algorithm  take  advantage  of  this  order.  In  this  note  we  shall  show  that. 
Indeed,  this  ordering  can  be  exploited  to  yield  an  algorithm  which  runs  In 
time  linear  In  the  number  of  the  edges.  Obviously,  since  each  edge  must  be 
examined,  the  time  of  our  algorithm  Is  optimal  within  a multiplicative 
constant . 

2.  It  Is  obvious  that  the  kernel  of  the  polygon  P,  being  the  Inter- 
section of  half -planes.  Is  a convex  polygon  K(P>.  We  shall  denote  P by  a 

doubly-linked  list  of  vertices  and  Intervening  edges  as  v.e>v,e,...v  ,e  .v,.. 

u u 1 1 n- i i u 

* This  work  was  supported  by  the  National  Science  Foundation  under  Grant 

NSF  MCS  76-17321  and  by  the  Joint  Services  Electronics  Program  under  Contract 
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We  also  impose  a direction  upon  each  edge  such  that  the  Interior  of  the 

polygon  lies  to  the  left  of  the  edge,  or,  equivalently,  the  boundary 

of  P is  directed  counterclockwise.  A vertex  v^  is  called  reflex  if  the 

angle  formed  by  its  two  adjacent  edges  e^  and  e^  meeting  at  v^  is  greater 

than  TT,  and  it  is  called  convex  otherwise. 

The  algorithm  we  shall  outline  scans  in  order  the  vertices  of  P and 

constructs  a sequence  of  polygonal  chains  K_ ,K, , . . . ,K  ,,  called  kernel  chains. 

0 1 n-1  “ 

Each  of  these  chains  is  a sequence  of  portions  of  straight  lines,  whose  first 
and  last  members  are  half-lines  and  all  others  are  line  segments.  As  we  shall 


show,  the  polygonal  chain  bounds  the  intersection  of  the  appropriate  half- 
planes determined  by  Cq ,e^, . . . ,e^.  Due  to  convexity,  the  angle  between  two 
consecutive  edges  of  a kernel  chain  is  always  < n.  Notationally , if  points 


i S£  i+1 


w and  w , belong  to  the  line  containing  the  edge  e of  P,  then  w « w . 
jL  1 I X S ^ X S ^ It  X 

denotes  the  segment  between  w,  and  w,  , and  directed  like  e ; moreover, 

i i+1  ®i 

A denotes  a point  at  infinity  and,  for  example,  Aew  denotes  a half  line 

terminating  at  vertex  w and  directed  like  edge  e. 

If  P has  no  reflex  vertex,  then  P is  convex  and  K(P)  “ P-  Thus 

let  Vq  be  a reflex  vertex  of  P.  Referring  to  figure  1,  we  set  Kq  equal 

to  the  intersection  of  the  half-planes  lying  to  the  left  of  edges  e^  j,  and 

Bq.  Notationally,  Kq  will  be  represented  by  the  string  of  symbols 

A e_  V,.  e , A*  For  each  K,  it  will  be  convenient  to  distinguish  two  vertices, 
U O n-1  1 

and  L^,  which  delimit  the  sequence  of  vertices  of  which  are  visible 
from  v^;  these  two  vertices  play,  as  we  shall  see,  a very  important  role 
in  the  construction  of  from  K^.  Obviously,  in  Kq  we  have  Fq  - Lq  ■ Vq. 
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Figure  1.  Illustration  of  kernel  chain 

We  now  develop  the  advancing  mechanism  of  the  algorithm,  i.e.,  the 
process  of  constructing  from  (K^,  F^,  L^) . For  later 

ease  of  reference,  It  Is  convenient  to  distinguish  a hierarchy  of  different 
cases. 

(1)  Vj^  Is  reflex  (see  figures  2a  and  b)  . In  this  case  lies  on 
or  to  the  left  of  the  half  line  v^^e^  and,  obviously,  Candidates 

for  are  only  points  belonging  to  the  subchain  delimited  by  F^  and  L^. 

We  now  examine  where  the  segment  lies  with  respect  to  ^e^ 

(1.1)  lies  to  the  right  of  Ae^  '^i+l  (^^8'**’®  2a).  We  scan  the 
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kernel  sequence  counterlockwise  from  F^,  until  we  find  a kernel  vertex  on 

or  to  the  left  of  '^i+1*  such  vertex  exists,  then  the  kernel 

is  empty.  Otherwise  we  can  determine  a unique  point  w'  as  the  intersection 

of  the  segment  w w and  Ae.  v and  set  F , «-  w'.  Next  we  scan  the 

kernel  sequence  clockwise  from  until  we  find  a point  w",  intersection 

of  Ae,  V..,  and  some  segment  w w Then  if  K.  * crw  ...w^  where 

i i+i  s s-1  1 s t“l 

a and  p are  sequences  of  alternating  edges  and  vertices,  we  set  = 

c^"e^w'p . 

(1.2)  lies  ot  the  left  of  Ae^  (figure  2b).  Let 


w » w . 1 

s s+1 


, . . . ,w^^^  be  the  sequence  of  the  kernel  vertices  between  F^  and 


L, , with  w = F,  and  w , = L. . Let  v.  denote  the  angle  measured  counter- 

clockwise  from  the  segment  (directed  from  Wj  to  to  e^^.  We 

successively  examine  the  angles  ^s+1’’’*’  '^*'til  we  find  some 

(0<  p S r),  such  that  is  minimal.  Notice  that  since  is  a 

convex  polygon,  only  w^,  > • • • >''s4.p»  '^g+p+i  need  to  be  examined  to 

find  w . . We  then  set  F,  w . and  K,  ♦-  K, . 
s+p  i+1  s+p  i+1  i 


(2)  v^  is  convex  (see  figures  3a,b,c,d).  In  this  case  F^^  lies 
on  or  to  the  left  of  the  half-line  v^^.  To  determine  ^i+1  , we 

distinguish  whether  the  vertex  L^  lies  to  the  left  of  v^e^^A  or  not. 

(2.1)  L^  lies  on  or  to  the  right  of  v^^e^A  (figures  3a, b).  We  scan 
the  kernel  sequence  clockwise  from  L^^  until  we  determine  a unique 
segment  w^w^-l  '^t  '^t-l  respectively,  to  the  right 

and  to  the  left  of  v^e^A:  we  can  then  determine  the  intersection  point 
w' of  distinguish  where 

to  w'.  Let  - Qfw^e. 


lies  with  respect 


(2.1.1)  w'e^A  (figure  3a).  Clearly  ♦-  v 

also,  we  obtain  - aw' 


(c) 


Figure  3,  Advancing  mechanism  when  v.  is  convex 
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(2.1.2)  3b).  Let  -y^  denote  the  counter- 
clockwise angle  from  the  directed  segment  w.v,.,  to  e..  If  w ...w  , 

® j 1+1  1 s s+r 

is  the  sequence  of  kernel  vertices  from  to  Lj^,  then  we  successively 

examine  the  angles  y , > until  we  find  a minimal  v , . We  then 

s s+l  s+p 

set  L.  . , ♦-  w'  , F w . , and  K.  , , - aw'e  ,A. 

i+l  1+1  s+p  1+1  1 

(2.2)  lies  to  the  left  of  v^e^A  (figures  3c, d).  Let  = 
aLj^e'A.  We  determine  the  intersection  w'  of  Le'A  and  v^^e^A- 

(2.2.1)  w'e^A  (figure  3c).  In  this  case,  we  set 


"i+l 


\+r  ^+1  S+l 


aL. e 'w ' e . V . . , e . A. 
1 1 1+1  1 


(2.2.2)  v^e^w'  (figure  3d).  In  this  case,  is 


determined  exactly  as  in  ttie  corresponding  case  described  in  (2.1.2) 

(figure  3b  ) whereas  ♦-  w'  and  — crL^e'w'e^A. 

In  all  of  the  above  cases,  it  is  immediate  to  realize  that  K 

the  intersection  of  K.  and  of  the  half-plane  to  the  left  of  e,. 

i i 


i+l 


is 


Using  the  advancing  mechanism  described  above,  we  ultimately 

obtain  the  kernel  chain  which,  if  K(P)  is  nonempty,  is  nonsimple 

(see  figure  4),  i.e.,  it  has  a crossing  point  w.  Our  remaining  task  is 

finding  w.  Let  K , = Ae'w.e!...  w e’A.  We  scan  the  edge  sequence  of 
n-1011  mm  w. 

starting  from  e^,  and  at  the  i-th  step,  for  i ::  2,  we  check  whether 

Wj^  lies  to  the  left  of  the  line  containing  e|,  directed  like  e!^.  Let  s 

be  the  smallest  value  of  i for  which  w^^  lies  to  the  right  of  e|.  Next 

we  scan  the  vertex  sequence  (wj^W2...)  until  we  reach  a vertex  w^,  such 

Chat  w , and  w lie  on  opposite  side  of  e'.  At  this  point  we  check 

whether  e'  and  e'  , intersect:  if  they  do,  their  intersection  is  the 
s r«l  ^ 

sought  w;  otherwise  we  replace  Wj^  with  w^  and  continue  the  process 


1. 

[ 


(repeating  the  alternate  scanning  of  the  edge  sequence  and  vertex 
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sequence)  until  the  intersection  is  found. 


Figure  4.  Finding  K(P)  from 

3.  We  now  analyze  the  performance  of  the  algorithm  outlined 

above . 

In  case  (1.1)  we  scan  starting  from  F^,  both  counterclockwise 
and  clockwise,  and  let  be  the  total  number  of  edges  visited  before 
finding  the  two  intersections  w'  and  w".  This  process  actually  removes 
- 2 edges  from  (those  comprised  between  w^  and  w^_j^  in  figure  2a) 
and  since  each  of  the  removed  edges  is  colinear  with  a distinct  edge  of 
P,  the  total  number  of  vertices  visited  by  the  algorithm  in  handling 
case  (1.1)  is  at  most  0(n). 

In  case  (1.2),  we  scan  counterclockwise  starting  from  F^^,  and 
clearly  (p+1)  is  the  total  number  of  vertices  visited  before  we  find 
But  in  this  process  the  distinguished  point  "F"  has  advanced  p positions 


( 


(from  F.  ■ w to  F. = w . ) counterclockwise.  Since  the  number  of 
vertices  of  any  is  at  most  0(n),  and  the  point  "F"  can  only  advance 
on  kernel  chains,  we  conclude  that  the  total  number  of  vertices  visited 
by  the  algorithm  in  handling  case  (1.2)  is  at  most  0(n). 

In  case  (2.1)  the  intersection  w'  of  v^e^A  and  involves 

scanning  clockwise  from  L^.  Let  be  the  total  number  of  edges 

visited  before  finding  w'.  This  process  actually  removes  edges  from 
(those  comprised  between  w^  and  A).  Here  again,  since  each  of  the 
removed  edges  is  colinear  with  a distinct  edge  of  P,  the  total  member 
of  vertices  visited  by  the  algorithm  in  finding  w'  in  case  (2.1)  is  at 
most  0(n). 

Case  (2.1.1)  requires  a constant  amount  of  work.  Case  (2.1.2) 
requires  globally  an  amount  of  work  at  most  0(n),  by  an  argument  identical 


to  that  developed  for  case  (1.2). 

The  discussion  of  cases  (2.2),  (2.2.1),  and  (2.2.2)  is  exactly 
analogous  to  that  of  (2,1),  <'2.1.1),  and  (2.1.2),  respectively. 

Finally,  it  is  straightforward  to  realize  that  finding  the  inter- 
section w in  requires  at  most  0(n)  operations. 

In  summary,  we  conclude  that  finding  the  kernel  of  a simple 
polygon  runs  in  time  0(n),  which  is  clearly  optimal  within  a factor. 
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